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Abstract. - We consider spatially separated qubits coupled to a thermal bosonic field that 
causes pure dephasing. Our focus is on the entanglement of two Bell states which for vanishing 
separation are known as robust and fragile entangled states. The reduced two-qubit dynamics 
is solved exactly and explicitly. Our results allow us to gain information about the robustness of 
two-qubit decoherence- free subspaces with respect to physical parameters such as temperature, 
qubit-bath coupling strength and spatial separation of the qubits. Moreover, we clarify the 
relation between single-qubit coherence and two-qubit entanglement and identify parameter 
regimes in which the terms robust and fragile are no longer appropriate. 



In recent years, we witnessed a fast development of the experimental realization of quan- 
tum bits and the maintenance of their quantum coherence. Despite this progress, decoherence 
stemming from a coupling of the qubits to a macroscopic environment remains a major ob- 
stacle for the implementation of a quantum computer. Several strategies are pursued to beat 
decoherence. One is dynamical decoupling: single qubits [1] or two-qubit gates [2] are effec- 
tively isolated from their environment by driving them with ac fields. Another strategy is 
quantum error correction, which requires a redundant encoding of a logical qubit by several 
physical qubits. Standard error correction protocols presuppose that all physical qubits couple 
to uncorrelated baths [3-5], which can be realized by putting qubits far apart [6]. A third, 
more direct strategy is the use of a decoherence- free subspace (DFS) . There, one logical qubit 
is encoded by several physical qubits, in such a way that the logical qubit states do not cou- 
ple to the environment [7-9]. Ideal DFSs occur when physical qubits couple via a collective 
coordinate to a common bath and are fairly robust against perturbations [10,11]. 

In view of the above, it is an important question whether different qubits are exposed to 
spatially correlated or uncorrelated noise. For example for charge qubits in quantum dots [12], 
a relevant source of decoherence is the coupling to substrate phonons [13,14]. Widely separated 
charge qubits at high temperatures will experience uncorrelated noise due to the phonon bath. 
The identification of a crossover regime at lower temperatures and smaller separations in which 
spatial correlations become important, requires a model in which the qubit separation and 
■ - - ■ the phonon bath are taken explicitly into account. 

In this Letter, we study the decoherence of spatially separated qubits which are coupled 
to a bosonic field that causes phase noise. On short time scales, phase noise is the main 
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source of decoherence in solid state environments [15]; on longer time scales, bit-flip noise 
becomes relevant as well [16], but we focus on the former herein. Phase noise is characterised 
by a qubit-bath coupling that commutes with the qubit Hamiltonian and, thus, allows an 
exact solution of the dissipative quantum dynamics. Still, the evaluation of the resulting 
exact expressions can be rather complex and often relies on approximations [3,6,7,17-22]. In 
contrast to those works, we here evaluate the exact solution in explicit form. We focus on the 
entanglement of two qubits prepared in particular Bell states and derive explicit expressions 
showing how their entanglement changes upon increasing their separation. For a wide range of 
parameters, we find that the dynamics is highly non-Markovian and that the entanglement can 
converge to relatively large values, even at high temperatures. We will discuss the robustness 
of the two-qubit DFS and clarify the relation between two-qubit entanglement and single-qubit 
coherence. 

Robust and fragile entangled states. Yu and Eberly [18] studied the entanglement 
dynamics for two qubits coupled to the same heat bath at identical positions, in particular 
for a preparation of the (maximally entangled) Bell states 

|01) + |10) |QQ) + |ii) m 

IV'robust) = \lp-) = ~q , ^fragile) = \lp+) = ~^ , (1) 

where |n) = \n l7 n 2 , ■ ■ ■ ,n N ) denotes the N qubit state with a vz \u) = (— l)™"|n), n v = 0,1 
and <j vz a Pauli matrix for qubit v. They found that the "robust state" lives in a DFS and 
consequently its entanglement is preserved, whereas the entanglement of the qubits initially 
in the "fragile state" decays to zero. In our study, we will also consider these initial states. 
The notation \ip±) has been introduced for writing equations more efficiently 

Two-qubit entanglement can be measured by the concurrence C[p] — max{0, ^f\\ — VA2 — 
v/As — \/Xa\ where denotes the eigenvalues of the matrix po-\y<JiyP* <J\y<Jiy in decreasing 
order; p* is the complex conjugate of p [23]. For maximally entangled states, one finds (7=1, 
while C vanishes for incoherent mixtures of product states. It was noticed [18] that for phase 
noise and two qubits at vanishing separation initially in either of the states , the concurrence 
is given by the absolute values of particular density matrix elements p min = (m|p|n), namely 
C_ = 2|p i,io| and C+ = 2|poo,ii|) respectively. These relations hold for spatially separated 
qubits as well. 

Qubits coupled to a bosonic field. - For modelling N identical qubits coupled to a homo- 
geneous bosonic environment, we employ a spin-boson Hamiltonian 

A N 

H = — <r«z + ^kblby + i? q _ b • (2) 

v=l k 

The first term denotes N qubits with energy splittings A and the second term represents a 
bosonic field with isotropic dispersion relation = u>k = c|k| and sound velocity c. The 
qubit-bath interaction 

N 

H q . h = hJ2^„ with & = ^ fffc (6 k e ik - x " +^e- ik - x "), (3) 

v=l k 

introduces phase noise due the coupling of qubit v via a vz to the field at position Xj,. We 
assumed identical and isotropic coupling strengths for each qubit, i.e. g^ v = gu- 

Furthermore, we assume that at initial time t = 0, the density matrix of qubit plus bath, 
i?(0), is of the Feynman- Vernon type, i.e. the bath is in thermal equilibrium and uncorrelated 
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with the qubits, i?(0) = p(0) ® p^ q , where oc exp{— ^ k fvjj\jD\.b^/ k^T} is the canonical 
ensemble of the bosons. Since we are exclusively interested in the behaviour of the qubits, we 
trace out the bath and obtain the qubits' reduced density operator p(t) = trb R(t). From the 
Liouville-von Neumann equation in the usual interaction picture, ihdR/dt — [i7 q _b(t), R], we 
obtain after some algebra for the density matrix elements the closed expression [6, 7, 20] 



p m ,„(i)=p m , n (0)e- A -"W +i ^ 



(4) 



The phases <^ m ,n(i) correspond to Lamb shifts which are brought about by time ordering. 
Since all quantities considered in the following are given by absolute values of single density 
matrix elements, these phases will not be relevant. The amplitude damping depends on the 
qubit separations ~x. vv > = x„ — x„' and reads 



A m ,„(i) 



cos(wfct) / hv k 
^ coth 



\2k B T 



N 

E 



[(-I)™* - (-l)"1e 



-ik-x,, 



(5) 



For the evaluation of these damping factors, it is convenient to introduce the spectral 
density j( d \uj) = ^ k gf<5(o; — cfc). For an homogeneous isotropic d-dimensional environment, 
it reads J^(uj) — aw(w/tj c ) <i_1 exp{— uj/uj c } [24], featuring the damping strength a and the 
cutoff frequency to c which for phonons is the Debye frequency. Then, the summation over the 
wave vectors k can be replaced by a frequency integration plus an integration over the solid 
angle. Evaluating the latter, we obtain for the concurrence the exact expression 

C { ±\t) = expj-8 duj J (d) M 

where G^(x) denotes a dimension-dependent geometrical factor which reads G^ x \x) = 
cos(x), G^{x) = Jq{x), and G^\x) = sin(a;)/a;, with J the zeroth-order Bessel function 
of the first kind. Note that G^ (0) = 1 for all dimensions, which causes the robust and the 
fragile behaviour of the Bell states Q for vanishing separation. For large argument, G^ 
decays for d — 2, 3, but not for d = 1. 

Below, we will compare the concurrence of a qubit pair with the single-qubit coherence 
|po,i(i)/poa(0)| = exp{— A^l(t)}, which we define for N = 1, i.e. when only one qubit is 
present. It is formally given by the rhs of eq. © but with the replacement 1 ± G^ — > |. 
Thus it is the geometrical factor G^ which determines the difference between single-qubit 
decoherence and entanglement decay of a qubit pair. In particular, the dimension dependence 
of C?w will turn out to be crucial but has been ignored in prior studies [20] . For qubits coupled 
to a three-dimensional bath, we find in the remote limit x±2 — > oo for both concurrences 
the relation C±(t) = exp{— 2A^l(t)}, which was also obtained for a model consisting of 
independent baths [21]. More generally, an intriguing corollary to eq. © is the exact relation 
C+\t) C^\t) = exp{— 4A^fl(t)}, which for arbitrary separation links the concurrences to the 
single-qubit coherence. It implies that if one of the concurrences vanishes, the single-qubit 
coherence must vanish as well. A finite single-qubit coherence, in turn, requires non- vanishing 
concurrences. 

In order to evaluate the concurrences © , we introduce the scaled time r = uj c t and transit 
time t 12 = lu c Xi2/c, and the scaled temperature 9 = k-QT/huj c . After inserting a Taylor 
expansion for cos(wr), we accomplish the frequency integrals in the resulting series [25]. Then 
we obtain for the single-qubit coherence and likewise for the concurrence an infinite product 



1 ± 



1 — cos(o;t) 



coth 



2k B T 



(6) 
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which can be combined into gamma functions and their derivatives. We restrict ourselves to 
the cases of one and three dimensions and find for the former case the single-qubit coherence 



r(0[l + ir])r(0[l-ir]) 



T 2 (0) 



4n 



1 



r 2 ^ 



2a 



and the concurrence 



c£>(r) 



-2A«(t) 



r(0[l-i(r 12 -T)])r(0[l-i(r 12 +T)]) 



r 2 (^[i-ir 12 ]) 



±8q 



(1 - in 2 ) 2 



±4a 



(7) 



(8) 



where T denotes the Euler gamma function. The corresponding expressions for a three- 
dimensional environment read 



-A< 3 J(r) 



exp | -4a ^20 2 Re|>i(0) - *i(0[l - it]) ] - 



t 2 (t 2 + 3) 

+ T 2 ) 2 



C£ 3) (r) = e" 2A o 3 i( T ) exp i±8a(J- Im[2* o (0[l - iru]) - *o(0[l - i(r - n 2 )]) 



-* o (0[l-i(r + r 12 )])] + 



t 2 (t 2 



'12 



3) 



(1 + t 2 2 )(t4-2t 2 [t 2 2 -1] + [1 + t 2 2 ] 2 ) 



(9) 



(10) 



where V^o & n d are Di-Gamma and Tri-Gamma functions, respectively. The importance 
of eqs. © and i|TT)|) lies in the fact that they explicitly yield the concurrences at all times 
for arbitrary spatial separations ri 2 , from a perfect DFS (ti 2 = 0) to uncorrelated noise 
(tl2 — * oo). In both expressions, the respective single-qubit coherences (JJJl and © appear. 

Time- evolution of the robust Bell state. - Let us first focus on the entanglement of a 
qubit pair that starts out in the robust state 1^-) and couples to a one-dimensional heat 
bath. Figure^ depicts the time-evolution of the concurrence for a temperature well below 
the Debye temperature. For vanishing separation, ti 2 = 0, the concurrence C^\t) remains 
at its initial value 1. This reflects the fact that then lives in a DFS and, consequently, is 
robust. For ti 2 > 0, we find that the concurrence initially decays until the transit time Ti 2 is 
reached. At time r = Ti 2 , the decay comes to a standstill and the concurrence remains at a 
finite value ci 1} (r -» oo) = (l + r 2 2 ) 4a |r(0[l - iTi 2 ])/r(0)| 16Q which becomes (1 + t 2 2 )- 4q for 
0^0. The time evolution allows the interpretation that before the transit time is reached, 
uncorrelated noise affects the qubits and entails an entanglement decay. After the transit time, 
the noise at the two positions is sufficiently correlated to establish a decoherence-free subspace. 
In the remote limit ri 2 — > oo, the concurrence of the robust state finally vanishes and, thus, 
the residual entanglement for finite ri 2 can be attributed to spatial bath correlations. We 
emphasise that for a one-dimensional bath, this scenario holds true for all temperatures. 

An intuitive physical picture for the observed entanglement dynamics is that at long times, 
decoherence is governed by the low-frequency modes of the bath. Owing to their large wave- 
lengths, these modes act effectively as a collective bath coordinate which leaves the entangle- 
ment robust. In higher dimensions, the role of the low-frequency modes is suppressed and, 
thus, the long-time behaviour may be significantly different [7]. 

Figure^ reveals that for a three-dimensional bath, in general the concurrence H10[) decays 
and saturates at a finite value which stays larger for closer qubits. In the experimentally 
relevant limit of low temperatures and ti 2 2 < 1, the final concurrence emerges as 
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Fig. 1 - Time evolution of the concurrence for the robust (left) and the fragile (right) Bell state for 
various spatial separations CT12/W. The qubits couple with a strength a — 0.01 to a one-dimensional 
(a,b) and a three-dimensional bath (c,d), respectively, at temperature 9 — kB.T/hu) c = 0.015. For 
qubit energies A = 0.01?kj c , the time range in the upper (lower) plots corresponds to 1.3 (0.03) 
coherent oscillations. 



However, the saturation generally occurs already at a time r <C T12, i.e. long before a field 
distortion can have propagated from one qubit to the other. To make this statement more 
quantitative, we numerically estimate the duration r* of the concurrence decay by the time 
at which 90% of the decay has happened. Figure shows that r* rs 1. In particular, t* 
is independent of the spatial separation T12, unless the qubits are very close. Hence the 
saturation of ci 3 ^(r) cannot be explained as a delayed build-up of a DFS. Instead, a single- 
qubit mechanism must be at work, since at times r < T12, the qubits experience effectively 
uncorrelated noise. This conjecture is supported by the resemblance of C_ (t) to the single- 

qubit coherence exp{— Aq \(t)} shown in fig. 01 A second difference to the one-dimensional 
case concerns the remote limit of the qubits: For T12 — ► 00, the stationary value is still finite. 
In this limit, in eq. © is negligible and the concurrence is given by the square of the 
finite single-qubit coherence. 

The concurrence (|llfl depends only quadratically on T12 and, thus, is rather robust against 
variations of the separation, provided the separation is small. Such robustness was pre- 
dicted [10] and confirmed experimentally [11] for symmetry-breaking perturbations. Interest- 
ingly, eq. {TJ| shows that the concurrence is robust against temperature variations as well, 
about which the theory in ref. [10] makes no predictions. 

Time- evolution of the fragile Bell state. - If the qubits are initially in the fragile state 
a one-dimensional bath causes an entanglement decay that becomes faster once the 
transit time is reached; see fig. As for the robust state cooperative effects only set 

in after a time 712. Whether the qubits are spatially separated or not, for a one-dimensional 
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Fig. 2 - Duration r* of the concurrence decay of the robust Bell state as a function of the 

transit time T12 for one- and three-dimensional baths. The dissipation strength is a = 0.01. 



Fig. 3 - Time evolution of the single-qubit coherence for various temperatures. The qubit is coupled 
with the strength a = 0.01 to a three-dimensional bath. 



environment their concurrence C + (r) ultimately decays to zero. The three-dimensional case 
again bears more surprises, as seen in fig.^i: the concurrence (|10f) decays, but the final value 

— > 00) is nonzero, in contrast to earlier statements [18, 19]. For low temperatures such 

— 1/2 

that 9 < 1, t 12 ' , the long-time limit is given by 

Cf (t - 00) = exp{-8a(l + ^ + ^) } . (12) 
This asymptotic value can be increased by reducing the temperature and by increasing the 

(3) (3) 

qubit separation. For a separation m > 1, C\ (00) = C_ (00), i.e. the concurrence of both 
the "robust" and the "fragile" state become identical; cf. the solid lines in figs.^^d. 

Conclusions. - The explicit evaluation of the exact reduced dynamics of two qubits with 
a nondemolitian coupling to a bosonic heat bath allowed us to investigate the consequences of 
a spatial qubit separation. We focused on two Bell states whose entanglement for vanishing 
separation is either robust or fragile. The most significant consequence of a finite spatial sep- 
aration is that the entanglement of the robust Bell state no longer remains robust: it decays 
initially, yet after a time t* — r* /u> c , it saturates. Indeed, it is interesting to find stable fi- 
nite bipartite entanglement even at high temperatures in our macroscopic system-bath model, 
which resembles recent results for a quite different model [26] . The duration t* of the entan- 
glement decay depends sensitively on the dimension of the environment: In one dimension, 
it equals the transit time of the field from one qubit to the other. In three dimensions, by 
contrast, the saturation is governed by a single-qubit effect and t* is approximately given by 
the inverse of the cutoff frequency u) c . Rather surprisingly, these durations exhibit only a 
weak temperature dependence. 

In a typical solid-state substrate, the Debye temperature is of the order 500 K which 
corresponds to Hlo c = 40 meV. For a qubit separation of 300 nm and a sound velocity c = 
3000 m/s, we find for a one-dimensional environment t* = t*/uj c « 10~ 10 s while in three 
dimensions, this time scale is much shorter: t* ~ 10~ 13 s. For a typical tunnel splitting A = 
10 fieV, the coherent oscillation period is 2ttH/A w 10~ 10 s, so that in the three-dimensional 
environment the concurrence of the robust Bell state undergoes a decay only during a very 
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short initial stage. At later times, a decoherence-free subspace is established and, thus, the 
concurrence stays robust. 

For the fragile Bell state, a three-dimensional environment in combination with a finite 
qubit separation prevents the entanglement from decaying entirely. If the qubits are sufficiently 
well separated, i.e. for X12 > c/co c , the entanglement of the "fragile" and the "robust" Bell 
state even assumes practically the same final value. In the above example, this is already 
the case if the qubit-qubit distance is larger than 1 fim, which usually holds for solid-state 
qubits. For typical parameters, the concurrence initially drops to and then remains at values of 
the order 0.9 already long before a first coherent oscillation is performed. Thus, uncorrelated 
phase noise creates decoherence-poor subspaces, which might be used for quantum information 
processing when complemented with quantum error correction protocols. 

* * * 
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